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The detection efficiency correction for the particle identification and kinematic fit selection 
criteria is investigeted for the branching fraction measurements in the V"' decays into wtt^tt", bnr, 



w/2(1270), ljK+K-, ujpp, 07r+7r-, <;/>/o(980) , (t,K+K- 



final states using 4 x 10 tp' data 



sample collected at BEPC/BES. Based on the corrected efficiencies, the preliminary results of these 
decay channels have been obtained, and the "12%" rule predicted by perturbative QCD theory 
tested. 
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Based on 4 x 10^ V' data sample collected at BEi- 
jing Spectrometer [||, BES collaboration measured 
the branching fractions of ip' decays into wtt+tt", 
&i7r, w/2(1270), ujK+K-, upp, ^tt+tt", 0/o(98O) 
, (j)K^K^, (j)pp final states. Perturbative QCD 
(PQCD) predicts ||^ that for any exclusive hadronic 
state h, the J/ip and ^p' decay branching fractions will 
scale as 



Qh = 



Biijy h) B{iJ;' 



B{J/iP h) B{J/tP 



12%, 



where the leptonic branching fractions are taken from 
the PDG tables This relation is known as the 
"12% rule". Although the rule works reasonably well 
for a number of specific decay modes, it fails severely 
in the case of the ip' two-body decays to the vector- 
pseudoscalar {VP) meson final states, pir and K*K 
This anomaly is commonly called the pir puzzle. 
In addition, the BES group has reported violations of 
the 12% rule for vector-tensor (VT) decay modes[||. 
Although a number of theoretical explanations have 
been proposed to explain this puzzle [Q , it seems that 
most of them do not provide a satisfactory solution. 
Therefore, it is meaningful to measure the branching 
fractions of these two- and three- body hadronic decay 
channels, and to test the " 12% rule" , for the improve- 
ment and development of existing models. 

The branching fraction of the decay process ip' ^ X 



* Translated from the paper which will be published in HEP & 
NP Vol.27(2),2003 (in Chinese) 
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X) = 



N4,,B{X ^Y)e{^P' ^ X 



where X stands for the intermediate state, Y the final 
state, e the detection efficiency, N^' the number of ip' 
events, and n"^" the observed signal events selected 
by applying certain creteria. The detection efficiency 
e is usually determined by Monte Carlo (MC) sim- 
ulation. Take the process of ip' ~^ loK^K^ as an 
example, the intermediate state X is ujK^K^ , the fi- 
nal state Y is tt^-k^ti^K^K^ . We produce events of 
ip' LoK^K^ TT^TT^Ti^K^K^ with phase space 
event generator HOWL, which then pass through the 
BES detector simulation package SOBER j|] and we 
acquire the MC simulated data. If we produce MC 
simulated data for iVg signal events, reconstruct them 
as same as for the real data and apply the same se- 
lection criteria for these reconstructed data, finally 
obtain N observed signal events, then, the detection 
efficiency is e^'^^=N/No. 

If the MC simulation and the reconstruction are suf- 
ficiently accurate, the detection efficiency determined 
in this way is correspondingly accurate. However, it 
is inevitable there exist much approximations in the 
simulation and reconstruction, therefore, the MC de- 
termined efficiency, e*^*-^, must differ from the actual 
value. If the difference is relatively small, it can be 
treated as a part of systematic error. If the devia- 
tion is relatively large, it is necessary to correct it to 
decrease the systematic uncertainty. 
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II. DEVIATION OF e^'^ AND CORRECTION 

The deviation between and actual efficiency 
depends on the event selection criteria. 

The event topology for the nine channels we investi- 
gated is either four prong (4P) or four prong plus two 

photons produced by a neutral pion decay (4P27). 
The general prc-sclcction criteria are following: 

1. The number of charged particles must be equal 
to four with not charge zero. 

2. The number of photon candidates must be equal 
to or greater than two for the decay channels 
containing 7r°. 

3. Particle identification (PID) 

For each charged track in an event, the Xpioi'^) 
and its corresponding Probpioii) values are cal- 
culated based on the measurements of dE/dx in 
the MDC and the time of flight in the TOF, with 
definitions 

xI/d(«) = X^E/dxii) + XTOpi}) 

Probpinii) = Prob{x]^^j^{i),ndfpiD), 

where ndfpjo = 2 is the number of degrees 
of freedom in the Xj^jp,{i) determination and 
Prohpinii) signifies the probability of this track 
having a particle i assignment. For final states 
containing pp, we require at least one of the 
charged tracks satisfy ProbpiD{p/p) > 0.01 > 
ProbpiD{TT/K), while for other channels ana- 
lyzed, the probability of a charged track for a 
candidate particle assignment is required to be 
greater than 0.01. For instance, in the selection 
of tj/ LoK^K~ TT^TT~TT^K~^K~ , wc rcquire 
that the probability of two tracks identified as 
pions must be greater than 0.01, and the prob- 
ability of another two tracks identified as kaons 
also be greater than 0.01. 

4. Kinematic fit 

A 4C (4 prong events) or 5C (4 prong plus two 

photon events) kinematic fit is performed for 
each event. To be selected for any candidate 
final state, the event probability given by the fit 
must be greater than 0.01. We will show that 
this requirement is uncorrclatcd to the PID re- 
quirement in our concrete event selection criteria 
(See Section V). 

5. The combined 



The combined x^, Xcorm defined as the sum 
of the values of the kinematic fit and those 
from each of the four particle identification as- 
signments: 

Xcom = X] Xp/£)(*) + Xfeine' 

which corresponds to the combined probability: 
Probcom = Prob{xlom^ndfcom), 

where ndfcom is the corresponding total num- 
ber of degrees of the freedom in the Xcom, de- 
termination. The final state with the largest 
Probcom is taken as the candidate assignment 
for each event. For example, in the selec- 
tion of V' ujK+K- Tr+TT-TT'^K+K- 
events, we require Pro6com(7r^7r~7r°if +if ~) 
must be greater than Probcom{'^~^'^~'^^'^~^''^~) 
and Probcom{K^K~'K^K^K~) to reject possi- 
ble backgrounds. 

The efficiency of PID requirement Probpioii) > 
0.01 is defined as ep/£)(z) = Ni/Nio, where Nio and Ni 
is the number of particle i before and after applying 
the PID requirement Probpioii) > 0.01. We wiU 
show in next sections that epfjj differs from fpfjj. 
Similary, the efficiency of kinematic fit requirement, 
e^^f differs from ef^, too. Therefore, we need to 
make adequate corrections for them. 

III. DEVIATION OF e^f^ AND 
CORRECTION 

If we select data samples of pure Tr,K,p,p par- 
ticles from real data by a set of criteria without 
PID requirement, these data samples should be un- 
biased to the PID requirement. By applying the cut 
Probpioii) > 0.01, we obtain epj^{i), which should 
be close to the real PID efficiency. The difference be- 
tween epjj-){i) and €pf^{i) reflects the deviation of the 
simulated MC efficiency from actual value. 

We select TT,K,p,p samples from processes f/)' — >■ 
7r^Tr~ J/ xjj, J/ ip Tr^Tr~Tr'^,K~^K~7r'^,pp7r'^,pp, be- 
cause they have large branching fractions and 4P and 
4P27 event topologies, which are as same as the pro- 
cesses we investigated. Therefore, the results should 
be applicable for our studies. 

• Determination of epjp,{n) 

epJj~i(Tr) is determined with charged pion sam- 
ple in tjj' TT^TT^ J/ip,J/-il) TT+TT^Tr" evciits, 
which are selected by following cuts: 
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1. Use cuts 1,2,4 stated in Section II to select 
ip' 7r+7r^7r'^7r+7r^ events 

2. Require Probfit{'K'^'iT^i:^iT^iT^) > 
Prob fit{'K^ Ti^ K^) to reject possible 
backgrounds 

3. Require \mfjj^i - 3.1| < 0.03 GeV to en- 
sure the existence of tt^tt^J/iP interme- 
diate state, where m^^J^^i is the invariant 
mass recoihng against the candidate tt'^tt^ 
pair 

Furthermore, we require one of daughter par- 
ticles from J/iJj decay selected with above cri- 
teria must satisfy Prohpio{'K) > 0.01, and 
ProbpiD{T^) > 10Probpi£i{K/p) to ensure it is 
a charged pion, then another track must also be 
a charged pion, with which we form a charged 
pion sample. This pion sample should be unbi- 
ased to pion's PID, because it is formed without 
PID requirement. By applying the requirement 
of ProbpioiT^) > 0.01 to this pion sample we 
can easily determine efficiency EpJoCtt). 

Fig. 1^ shows the dE/dX distributions of candi- 
date charged pions selected with above criteria 
for MC and real data. In the Figure, the pi- 
ons with higher momentum are produced from 
J/ip decays, while those with lower momentum 
are the "direct" pions from i/;' tt+tt" J/^/; de- 
cays. We see that the purity of pion sample is 
good, no pollution from K^p,p seen. Also, the 
consistency between MC and real data is nice. 




P(Gev/c) P(Gev/c) 

FIG. 1: dE/dX vs. p (momentum) of charged pions in 
-K^ -K^ J / ip ,J / ^jj —* TT^TT^Tr" events 

• Determination of epf]^{K) 

epjp, {K) is determined with charged kaon sam- 
ple in -0' — » 7r+7r~ J/'i/', J/V' — > K^K~-k^ events, 
which are selected by following cuts: 



1. Use cuts 1,2,4 stated in Section 2 to select 

— > Ti^TT^Ti^K^K^ events 

2. Require Prob fit{Tr^Tr^Tr^Tr^Tr^) < 
ProbfitiTT'^TT^Tr'^K^K^) to reject possible 
backgrounds 

3. Require \m'^^J'^'^i - 3.1| < 0.03 GeV to en- 
sure the existence of n^n^ J/ijj intermedi- 
ate state. 

In addition, we require one of daughter particles 
from J ftp decay must satisfy Probpip){K) > 
0.01, and ProbpiD{K) > lOProbpiDin/p) to 
ensure it is a charged kaon, then another track 
must also be a charged kaon, with which we 
form a charged kaon sample. This kaon sample 
should be unbiased to kaon's FID, because it is 
formed without FID requirement. By applying 
the requirement of Probpjp){K) > 0.01 to this 
kaon sample we can easily determine efhciency 

• Determination of epj^{p/p) 

^pfoip/p) determined with p/p sample in 
ip' — > TT~^TT~ J/ip,J/ip — > ppn'^ events, which are 
selected by following cuts: 

1. Use cuts 1,2,4 stated in Section II to select 
tp' 7r+7r^7r°pp events 

2. Require Prob fit{Tr^TT^ n'^pp) must be 
greater then Pro6/if (7r+7r^7r'^7r+7r^) and 
Probfit{TT~^TT~TT'^K~^K^) to reject possible 
backgrounds 

3. Require Im'^^^^ - 3.1| < 0.03 GeV to en- 
sure the existence of vr+vr^ J/i/j intermedi- 
ate state. 

We further require one of daughter particles 
from J/ip decay must satisfy Probpjr){p/p) > 
0.01, and Probpioip/p) > WProbpiDiTr/ K) to 
ensure it is a proton (antiproton) , then an- 
other track must be an antiproton (proton), 
with which we form p/p sample. This sample 
should be unbiased to {p/p)''s FID, because it 
is formed without FID requirement. By apply- 
ing the requirement of Probpip)(j)/p) > 0.01 to 
this sample we can easily determine efficiency 
^'^Jnip/P)- 

Fig. ^ shows the dE/dX distributions of p/p sam- 
ples selected with above criteria for MC and real data. 
In the Figure we see that the purity of proton sample 
is good, the consistency between MC and real data is 
nice; while the antiproton sample is polluted by other 
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particles, which are produced by the annihilation of p 

in detector material and not correctly simulated in the TABLE I. epiD(7r),epj£,(7r) and /piD(7r) 

Monte Carlo, therefore the MC determined efficiency 
has to be corrected. With these p/p samples we deter- 



p(GcV/c) 6gfo(^) e^fnin) fpinjTv) 

_ < 0.5 0.950 ± 0.004 0.966 ± 0.002 0.983 ± 0.005 

mine ep^^(p/p), and in the determination of ep^^(p) o.5 -> 0.8 0.979 ± 0.022 1.000 ±0.001 0.979 ± 0.021 
only the p sample which dE/ dX values located on the g.S ^ 1.2 0.961 ± 0.014 0.984 ± 0.009 0.977 ± 0.015 
expected band is selected to cut out other particles 1.2 -> 1.4 0.956 ± 0.017 0.989 ± 0.005 0.967 ± 0.017 
pollution. ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 



TABLE IL e^iD,e'¥fo and fno for K,p,p 




P(Gev/c) 



0.5 1 1.5 

P(Gev/c) 



FIG. 2: dE/dX vs. p (momentum) of p/p 
Tv'^Tv~ J/ ip, J/ ip ppty" events 



in ^' 



Table I,II list the e|?f^(i) and e'),%{i) for PID re- 
quirement ProbpiD > 0.01, and the corresponding 
correction factor fpioii) = ^pJoi^) / ^pFoi''')- From 
tables we know that, fpiD(T^) and fpioiK) are close 
to 1, while fpioip) and especially fpioip) deviate 
from 1 significantly. Therefore, for the final states 
containing pp pair we will not require that both 
Probpioip) > 0.01 and ProbpiD{p) > 0.01 must be 
satisfied, instead, we just require one of them be sat- 
isfied. In this case the correction factor fpioip/p) is 
close to 1, as shown in the last row of Table II. Be- 
cause p and p has to be produced in a pair, if one of 
them is correctly identified, its twin is also correctly 
identified. 

Take tp' ipK+K- K+K~K+K- as an exam- 
ple, if we use PID requirement Probpi£i{K) > 0.01 
to select four candidate kaons, the correction factor is 



Particle i 






fpiD {i) 




0.952 ± 0.020 


0.980 ± 0.004 


0.971 ± 0.020 


P 


0.889 ± 0.035 


0.936 ± 0.006 


0.950 ± 0.030 


P 


0.759 ± 0.052 


0.929 ± 0.006 


0.817 ±0.043 


p/p* 


0.973 ±0.016 


0.996 ± 0.008 


0.977 ±0.018 



f PI Dim = IfPIDiK)]' 



* either p or p satisfies PID requirement 



IV. DEVIATION OF e 



fit 



AND CORRECTION 



particle reaction should conserve 4-momcntum in 
initial and final states. In experiment, the measured 
4-momentum of particle will somewhat deviate from 
its true values due to measurement errors. Therefore, 
the 4-momentum conservation is not accurate for mea- 
sured values. The kinematic fit is an algorithm, which 
uses least square priciple, by varying measured 4- 
momenta of particles to satisfy the 4-momentum con- 
servation and other physics constraint (for instance, 
some particles in final state coming from an inter- 
mediate resonance decay) to the possibly best level 
within measurement errors. Kinematic fit is a power- 
ful tool for signal events selection and background re- 
jection. In our studies we use requirement Probfu > a 
(0 < a < 1) to select signal events. However, because 
MC could not simulate physics reality without error, 
this requirement will cause bias of ^Yi^ ■ 

If we select a data sample of certain events with- 
out kinematic fit, this sample should be bias free to 
the kinematic fit requirement. By applying Probfu > 
0.01 to this data sample, we acquire ef^^, which should 
be close to the true efficiency. Then, the difference 
between and ejl^ refiects the deviation of MC 
simulation from real data due to the kinematic fit re- 
quirement. 

We select ip' tt+tt" J/-0, J/ip pp events with- 
out using kinematic fit requirement by following cuts: 
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1. < 0.6 GeV/c (in this process, the momem- 
tum of pions are rather low). 



2. Require \m 



'recoil 



'3.1\< 0.03 GeV to ensure the 
existence of tt^tt^ J/ip intermediate state. 

3. require 0.9 < Pp/p < 1.5GeV/c to reject possible 
J/ip e+e^,//+/i^ backgrounds. 

4. Require Prohpioip) > 0.01, ProbpiD{p) > 
0.01, to ensure ip' tt^tt^ J/ip,J/Tp pp 
events selected. 

Based on the selected events we determine the branch- 
ing fraction oi J/ip ^ pp to be 2.12 x 10~^, consis- 
tent with the PDG value Q , which indicates the reli- 
ability of events sample. Applying the requirement 
Prohfit{-K^Tr'pp) > 0.01 obtains ef^{TT^n^pp) ~ 
0.773 ±0.015, while from MC we get €jlf{TT+n-pp) = 
0.948 ± 0.08, therefore, the correction factor is 
ffu{-n+TT-pp) = 0.815 ± 0.017. 

Reference has shown that the correction fac- 
tor for kinematic requirement Probfit(TT^'!T~ fj,'^ fj,^) > 
0.01 is ffiti'n-'^TT'fi'^^^) = 0.859 in process V' ^ 
■n^ t:^ J / %p ,J / ij} /^^M^. 

Notice that, in the momentum region involved 
under our study, muon has the smallest, and pp 
the biggest interaction with the detector material in 
/i, TT, K,pp four type particles, therefore, we use same 
correction factor //jt(4P) = 0.85±0.05 for all 4F final 
states in our study. 

e^e^ — > 77 events have distinctive topology (two 
oppositely directed neutral tracks with the energy 
close to the beam energy, etc.), hence it is easy to 
identify and select with no need of kinematic fit re- 
quirement. After the correction to the error matrix of 
neutral track it is found //it(e+e~ — > 77) = 0.996 
for the kinematic fit requirement Probfit{e'^e~ — > 
77) > 0.01. This indicates that the existence of neu- 
tral tracks does not affect the correction factor of kine- 
matic fit requirement. Therefore, to the 4^27 final 
states under study, we use the same correction factor 
as to the 4P final states for kinematic fit requirement, 
but the uncertainty is increased from 0.05 to 0.08, 
namely, //,t(4F27) = 0.85 ± 0.08. 



as the function of a for kinematic fit requirement 
Probfit > a, in which the point at a < denotes 
the kinematic fit is not applied. The statistical error 
is expressed by error bar. From the Figure we see 
that ebias (therefore fpiniT^)) is basically a constant 
for different a within errors. This implies that the cor- 
rection of PID is basically irrelavent to the correction 
of kinematic fit in our case, the correlation between 
them is rather weak, hence can be negligible. 



0.02 0.04 0.06 0.08 0.1 



FIG. 3: Correlation between e*ff (tt) and epfij{n). See 
text for the meaning of ei,ias and a. 



VI. SUMMARY 

Based on the method described above we determine 
the efficiency correction factors for nine final states 
listed in Table III, with which we give the prelimi- 
nary results of their branching fractions and corre- 
sponding Qh values (the ratio of branching fractions 
of -0' and J/ip) ijl^. Among these, the branching 
fractions for ti;/2(1270) and bfn^ supersede previous 
BES results |ll| ; while all other branching fractions 
are first measurements for these decays. The ratio of 
ijj/2(1270) (VT mode) is suppressed by a factor of five 
with respect to the PQCD expectation, and those of 
clitt+tt", Lopp, and (pK'^K" are suppressed by a factor 
of two, while those of 6f tt^ (AP mode), 0/o(98O) (VS 
mode), wif+if" and (/)7r7r are consistent with PQCD 
expectation within errors. As to the (ppp channel, we 
need more statistics. 



V. CORRELATION BETWEEN efto AND ef^ 

With the pion sample selected from ?/;' 
tt+tt^J/t/;, J/ip fJ-^ 1^^ events we determine 

epJo(7r), and using the MC simulation for same pro- 
cess we determine ep^/f,(7r). Difine e^as = '^Mci^l — 
1 = fpioiT^) - 1- Fig- a shows the euas value 
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TABLE III: Efficiency correction factors 



Channel 



Effi. Corr. Factor (%) 



cj/2(1270) 

ujK+K- 

uipp 

0/o(98O)(/o -*7r+7r- 
(j)pp 



77.4 ±8.7 
77.4 ± 8.7 

77.4 ± 8.7 
78.0 ±8.3 
85.0 ±8.0 
75.6 ±7.6 
75.6 ±7.6 

76.5 ±7.7 
85.0 ±5.0 
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